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The relation between angular diameter distance and redshift {Ja-z relation) in a spherically 
symmetric dust-shell universe is studied. This model has large inhomogeneities of matter distri¬ 
bution on small scales. We have discovered that the relation agrees with that of an appropriate 
Friedmann-Lemaitre(FL) model if we set a “homogeneous” expansion law and a “homogeneous” 
averaged density held. This will support the averaging hypothesis that a universe looks like a FL 
model in spite of small-scale fluctuations of density held, if its averaged density held is homogeneous 
on large scales. 


I. INTRODUCTION 

The standard big bang model is based on the assumption of homogeneous and isotropic distribution of matter 
and radiation. This assumption then leads to the Robertson-Walker (RW) space-time geometry and the Friedmann- 
Lemaitre (FL) universe model^ through the Einstein equations. This model has succeeded in explaining various 
important observational facts: the Hubble’s expansion law, the contents of light elements and the isotropy of cosmic 
microwave background radiation (CMBR) ||]. 

The CMBR conversely gives a strong observational basis for the assumption of homogeneity and isotropy of our 
universe by its highly isotropic distribution together with the Copernican principle. Indeed, the deviation of our 
universe from a homogeneous and isotropic space is as small as 10“® at the stage of decoupling. Thus our universe 
is well approximated by a FL model before this stage. On the other hand, the present universe is highly inhomogeneous 
on small scales; the density contrast is of the order of 10^° in the solar system, 10® on galactic scales, and of the order 
of unity even on the scale of superclusters. We have to go beyond FL models in considering such systems. 

For a long time we have regarded that a FL model is a large-scale “average” of a locally inhomogeneous universe 
(averaging hypothesis). Even though the observational data are consistent with the picture that our universe is 
described well by a RW metric, we are not sure how to derive the background FL model from the inhomogeneous 
universe by any averaging procedure, or how the non-linear inhomogeneities on small scales affect large-scale behavior 
of the universe |^. Although one can derive a background FL model from observations of the nearby galaxies with 
any rule of averaging one likes, it is uncertain whether this background FL model agrees with the FL model whose 
cosmological parameters was defined at the stage of decoupling. The discrepancies might appear, for example, in the 
estimate of the density of baryonic matter, the density parameter, the age of our universe, and so on. These still 
remain a non-trivial question to which we have to give a clear answer. 

A number of approaches have been made to study how the small-scale inhomogeneities affect the global dynamics 
when averaged on larger scales ill- The first work which explicitly showed the existence of such an effect was 
performed by Futamase Q|. Assuming small deviation of the space-time geometry from the RW one, Futamase 
constructed an elegant formalism by the use of the Post-Newtonian expansion and Isaacson’s prescription to take 
into account the back reaction of the small-scale inhomogeneities on the global cosmic expansion. After his works, 
Buchert and Ehlers studied this back reaction problem in the framework of the Newtonian cosmology in which 
the corresponding background FL model is uniquely determined through the spatial averaging of physical quantities 
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without any uncontrolled approximation. They showed that inhomogeneities do not influence the overall expansion 
in spatially compact models (the topology of its spatial section is T^) if they are averaged over the whole space. 
In other cases, the significance of inhomogeneities may depend on the cosmology one adopts. Another interesting 
approach has been made by Carfora and Piotrkowska |Q, in which 3-dimensional geometry is deformed according to 
the Ricci-Hamiltonian flow which they say would be equivalent to changing the scale of averaging. They derive a 
homogeneous geometry which corresponds to a large-scale average of an inhomogeneous universe, and discussed the 
effect of inhomogeneities on the cosmological parameters. 

In spite of these works, the effect of inhomogeneities on the cosmological parameters remains unclear; there are 
even apparent discrepancies in their statements. The reason seems to lie in the different definition of back reaction of 
inhomogeneities 0; its definition is ambiguous since they do not treat observable quantities in these works. Thus, 
in order to understand clearly the effects of inhomogeneities, it is necessary to relate them with physical quantities 
which we can give an unambiguous definition. 

Bildhauer Q showed that the global cosmic expansion rate is not isotropic if the back reaction of small-scale 
inhomogeneities are taken into account and then investigated its observational effects on the CMBR. After this work, 
Bildhauer and Futamase (D discussed the possibility that the observed dipole anisotropy in the CMBR comes from this 
effect. These works are significant in the sense that the back reactions on the observable quantities were discussed, but 
they did not consider the cases where we try to determine the cosmological parameters by observing an inhomogeneous 
universe. 

In this paper, we investigate the relation between distance and redshift in an inhomogeneous universe. As an 
inhomogeneous model, we take a spherically symmetric dust-shell model, in which light rays travel from each dust 
shell toward an observer at the center. We calculate angular diameter distance-redshift (cIa-z) relation and compare 
it with that of a FL model. Our main goal in this paper is to clarify the condition under which the distance-redshift 
relation of a dust-shell universe behaves like that of a FL model. We investigate the behavior of the averaged density 
around the observer at the center when we gradually extend the averaged region to the outer shells, and discuss 
the relation among the behavior of the averaged density, the energy density of the FL model, and the observed 
distance-redshift relation. 

We note two interesting points of a dust-shell model. First, its dynamics is exactly solved; it is not necessary to 
assume the existence of homogeneous background in order to obtain the behavior of density fluctuations. Secondly, 
it can treat a discrete mass distribution where the linear perturbation theory is invalid, and can also treat a highly 
general-relativistic situation where the scale of inhomogeneities are comparable to the horizon scale. 

The organization of this paper is as follows. In the next section, we present basic equations for the dynamics of a 
dust-shell universe, distance to the shells from the center, and redshift of the shells measured by an observer at the 
center. We give our results and discussion on cIa-z relation and averaged density in section 3, followed by concluding 
remarks in section 4. 

We follow the sign convention of the Riemann tensor and the metric tensor in and adopt the unit of c = 1. 


II. FORMULATION OF DUST-SHELL UNIVERSE 


A. Equation of motion of dust—shell 


We put a number of spherically symmetric shells whose centers are common at r = 0 (see Figure 1). The innermost 
shell is called the first shell, the next one is called the second shell, and so on. The region enclosed by the (z — l)-th 
shell and z-th shell is called the z-th region. Each shell is infinitesimally thin, characterized by the surface stress-energy 
tensor which is given by 

5“'' = lim J dx , (2.1) 


where x is a Gaussian coordinate (x = 0 on the shell) in the direction normal to the shell. 

Since each region between the shells is vacuum, the space-time is described by the Schwarzschild geometry by the 
Birkhoff’s theorem. The line element in the z-th region is written in the form 


dsf 



dt^ -|- 



2Gmi^ 


dr^ +r^dn^, 


( 2 . 2 ) 


where the parameter mt will be referred to as a gravitational mass (mi = 0), and dfl is the line element of a unit 
sphere. 
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We first derive the expansion law of the i-th shell following Let n“ be a unit space-like vector normal to 

the trajectory of the shell, and define the projection operator — UaU^■ From the projected Einstein equation 


RabKK = SttG ( Tab - ^gabT ) 


(2.3) 


one obtains 


£nKad + ^Rcd - KKad = 8^G ( TabKh'', - ^h^dT 


1 . 


(2.4) 


where £n is the Lie derivative along n“ and ^Rcd is the 3-dimensional Ricci tensor on the timelike hypersurface 
generated by the motion of the shell. The extrinsic curvature Kab is defined by Kab = —^h^hf£nhcd, and K = K^, 
T = Ta - Integration of equation (2.4) over an infinitesimal range along yields 


~ ^nh = SttG f Sab - -habS 


^ab 


(2.5) 


wher e the suffix ‘-I-’ denotes a quantity evaluated at the outside of the shell and ’ at the inside. Using equation 
(2.5) and the Gauss-Codazzi relation 2Gab'n£n^ = —^R + KabK^"^ — K^, one finds that the following relation holds 
for a dust-shell: 


S‘^'{Kb + Kab)=Q ■ 


( 2 . 6 ) 


obtain the following equation for the circumferential radius (the “expansion law” of the dust shell): 


Combining equations (2.5) and ( |2.6D and substituting the expression of the metric and the stress-energy tensor, we 

-0 + " , (2.7) 


drj 

dr 


2GM+{i) 


M{i) 


s(*) 


4rf 


where M± is defined by 


M+(i) s 


( 2 . 8 ) 


M_{i) = rrii+i - rrii , (2.9) 

and T is the proper time of the shell. We have also introduced the “baryonic” mass ms{i) given by ms{i) = Ansirf 
where Si is the surface density of the i-th shell; Si = —Sa{i)- It can be shown this baryonic mass is constant of 
motion by the conservation law, S'^.^ = 0 where the semicolon denotes the 3-dimensional covariant derivative on the 
trajectory of the shell. 

In this paper, we shall use a common proper time r for all the shells. The relation between r and the time 
coordinate, t, in Schwarzschild space-time is obtained as follows. First, note that two Schwarzschild time coordinates 
are assigned to each shell: the i-th shell has the time t(-)i measured in the f-th region and t(+)i measured in the 
(i -|- l)-th region. Then, from the normalization condition for the 4-velocity of a shell, we obtain 

( 2 . 10 ) 


( 2 . 11 ) 

The procedure to determine the initial points of each time coordinate is described later. 


dt 


i+)i 


dr 

dt( — )i 
dr 


n - 2Gmi+i 


Ti — 2Gmi 


1 - 


2Gm. 


■i+l 


drj 

dr 


1 - 


2Gm, 


n 


dvi 

dr 
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B. The solution and initial condition 


In order to specify a dus t-sh ell universe, we have to fix the parameters in equation (2.7) and the initial hypersurface. 
We first rewrite equation (^) in the form corresponding to the Hubble equation of FL models. We denote the initial 
circumferential radius of the z-th shell by cci, i.e., 


on initial hypersurface. 


( 2 . 12 ) 


We define pi by 




(2.13) 


and ki by 


M_(i) ^ ^ 
m^ii) 


= 1 — kiX~ 


Then the expansion law of the dust-shell is written as 


(2.14) 


-p-) 

n dr J 3 \r,J Xn 


+ 


(z) 

4rf 


(2.15) 


We see that the first term behaves like a non-relativistic matter term in the Hubble equation of FL models, the 
second and the third like a curvature and a radiation source term. From this point of view, pi and ki play roles of 
the “energy density” and the “curvature”, respectively. The radiation-like term might be regarded as the effect of 
the binding energy of the shell |Q . Further it is worthwhile to note that this radiation-like term is consistent with 
Futamase’s result about the effect of the small-scale inhomogeneities on the global cosmic expansion |^. Seeing this, 
one may expect that the inhomogeneities tend to make the Hubble parameter larger compared with a homogeneous 
universe which has the same “energy density” of non-relativistic matter. However, this radiation-like term does not 
necessarily imply the larger Hubble parameter. In order to see the effect of this term on the Hubble parameter, we 
need to investigate the distance-redshift relation by solving the null geodesic equations and compare the result in the 
inhomogeneous space-time and that of the FL model. Such an analysis will be per form ed in the following sections. 

A dust-shell universe is specified if we set the parameters contained in equation ( 2.15 ), i.e., pi, ki, Xi, and an initial 
hypersurface. When we increase the number of the shells to infinity with pi and ki being finite and independent 
of z (we will mention this limit as “large N limit”), the dust-shell universe approaches a FL model if we take an 
appropriate initial hypersurface, as we will see in section 3. Then the parameters pi and ki agree with the ordinary 
energy density and curvature in the Hubble equation. 

We take 


ki = 0 


(2.16) 


for all z in the remaining sections. This means that the increase in the gravitational mass rrii is equal to the baryonic 
mass of the shell. In other words, the kinetic energy of the shell balances with the potential energy, and hence the 
total energy becomes equal to the rest mass. This is also the simplest case which approaches a flat FL model in the 
large N limit. Choice of the other parameters and the initial hypersurface will be discussed in section 3. 

For notational convenience, we shall introduce the following quantities: 


Pi = 2Gmi, Vi = -Gmsii) and Ui = 2GM+(z) 


From equations (p.q) and (^), the following relations hold 


Pi — Oi 2vi , 

Pi-\-l — O'i 4 - 2vi. 


(2.17) 


(2.18) 

(2.19) 


Then the equation for the circumferential radius, of the z-th shell is written in the form 


dvi 

dr 


I % 


( 2 . 20 ) 
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where we have assumed that each shell initially expands. From equations (2.10) and (2.11), the equations for the 
Schwarzschild time coordinates, t(±)i, are obtained as 


dt 


(±)* 


n T Vi 


dr 




Thus, from equations (2.20) and (2.21), the equations for the relations between t(±)i and are given by 

Tiivi Ui) 

dn (n - aiT ‘2vi)y/(Jiri + vf 

The above equations can be integrated easily to give the solution for > Ui ± 2vi in the form 
^(±)i(^i) = ± 2^'i) In —--^-- + r(±)i(rj) + T(±)i, 


+ (cTj ± Vi) 


where 


T{±)i{T’i) = + 3cr,^ ± Sr-iCTi - 2i/f) aiTi + ; 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


and 7(j-)j’s are integration constants. 

There is a coordinate singularity on the Killing horizon; t(±)i becomes infinite on n = Ui ± 2vi. For further 
calculation the null coordinate is convenient since we are interested in the null geodesics in this space-time. Hence, we 
shall adopt the Kruscal null coordinates which has no coordinate singularity. Outside the horizon in the i-th region, 
r > fii, the Kruscal null coordinate is given by 


U = -2v^(r - Mi)" exp . 

V = +2y^(r - fii)^ exp > 


(2.25) 

(2.26) 


where U and V correspond to the retarded time and the advanced time. Using these Kruscal coordinates, the line 
element in the i-th region is expressed as 


dsi = — ^exp ( — — ) dUdV 


2r 


fJ-i 


^dn\ 


(2.27) 


Similarly to the Schwarzschild time coordinate, two pairs of Kruscal null coordinates, and V(±)i, are assigned 
to each shell. Using equation (2.23), we obtain the Kruscal null coordinates labeling the i-th shell in the form 


U(^±),{n) = -2V2Jl± 


2vi 


Ui 


cTiTi + vf + ai±Vi\ exp 


^i - '^(±)i - T'{±)i 




2vi 


(±)i(u) = +2^2^/ \±-^[ \lGin + vf - CTi T O ) exp 


2(cri ± 2vi) 

+ ’^(±)i + ^(±)i 


2(cri ± 2vi) 


(2.28) 

(2.29) 


As expected, and V(^±)i are finite on = (7^ ± 2vi and are well defined also for < cji ± 2vi. When Ui ± 2vi is 
larger than r^, both and V(±)i are negative. This means that the i-th shell with < Ui ± 2vi is located in the 

white hole part of the Schwarzschild space-time. This situation occurs for the shells beyond the horizon scale. 

The determination of the integration constants T(^±)i corresponds to the choice of the initial hypersurface. The 
procedure to construct the initial hypersurface we adopt is summarized as follows; First, we choose a unit spacelike 
vector which is directed outward in the ordinary sense with respect to r. Taking this vector as a starting tangential 
vector, we extend a spacelike geodesic curve until it reaches the second shell. This spacelike geodesic curve defines 
the simultaneous hypersurface in the region between the first shell and the second shell. Next we extend from this 
intersection towards the third shell another spacelike geodesic which starts from another spacelike vector at the second 
shell. This second spacelike geodesic generates a spacelike hypersurface in this region. Repeating this process, we 
complete the whole initial hypersurface. 

From the above procedure, the integration constant of equation (2.23) is determined as follows. In the z-th region, 
the tangent vector of the spacelike geodesic is written as 
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l + E^ 


IM 

r 


(2.30) 



and R=x 

V r J 

V 


and the other components vanish, where Ei is an integration constant associated with the geodesic equation and will 
be determined by the condition which we will see in the next section. From equation ( 2.30| ), the equation for the 
trajectory of the spacelike geodesic in the {t, r) plane is given by 


dt Eir^ 

{r - ti^)^/{l + Ef)r- n. 

Integrating the above equation, we obtain 


t = Fi{r) + Di, 


where 

1 + E" \s/(l + Ef)r-n. + E,^) 

Ei{3 + 2El)fj,, ln(y(r+ Ef)r - + ^(1 + Ef)r) 

{1 + Efr/^ 


(2.31) 


(2.32) 


(2.33) 


and Di is an integration constant. Initially, we set t{+)i = t{-)i and t(±)i = 0. Then, since = Fi-i{xi-i) + Di-i 

and t(+)i-i = Fi{x^-l) + A, we find 

Di = A-i + F^_i{xi-i) - Fi{xi-i), (2.34) 

and Di = —Fi(a^i). From the above recurrent relation, we obtain 

i-l 

D^ = -F,{x,-i) + [Fj{xj) - Fj{xj-i)] . (2.35) 

t=2 

Using these relations, we obtain the integration constants, for i > 2 as 


^±)i ~ 

i=2 


/ , o M ± Iy^) 

- (cr, ± 2iyi) In -^^-- - T(^±)i{xi), 

^UiXi + vf + (di ± Ui) 


and for i = 1 as 


(2.36) 


7(±)i = -(cti ± 2vi) In 


yJdlXi + - (gi ± Vl) 

yjuixi + vl + (tTl ± Vl) 


T(±)i{xi) 


(2.37) 


C. Redshift and diameter distance 

We consider a light ray which is emitted from each shell toward an observer at the center. The light ray goes along 
a future directed ingoing radial null geodesic, where “ingoing” refers to the direction from a shell toward shells labeled 
by a smaller number. 

An ingoing radial null geodesic is specified by a constant V in the Kruscal null coordinate. The circumferential 
radius of the Fth shell when it intersects the null geodesic is denoted by Ri. The outermost shell considered here is 
labeled by M. Then Rm = xm and hence in the M-th region, V = V(-)m(xm) is satisfied along the null geodesic. 
Thus, on the {M — l)-th shell, the following relation holds 

y(+)M-i{RM-i) = V(_)m{xm)- (2.38) 

This equation determines Rm-i- We obtain the circumferential radii of all the shells at the intersection with the null 
geodesic by the same procedure, i.e., by solving the following recurrent relation: 
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y{+)i{Ri) — V(-)i+l{Ri+l)- 


( 2 . 39 ) 


We can determine Ri from the given through this equation. 

In order to derive the expression of redshift, we first write down the components of the null geodesic tangent in the 
j-th region, k^(i), which is given in the Kruscal null coordinate as 

(i) = — exp ( — j uji, (2.40) 

Mi Kk-i / 

and the other components vanish, where Ui is an integration constant associated with the geodesic equation. We 
require that the observed frequency of the photon at each shell is uniquely determined. The observed frequency, 
LUohii), at the i-th shell is given by 


UJoh{^) = 


uJoh{i + 1 ) = -k^{i)u'^_^{i + 1 ) = -uJi 


1 '^l^(-)i+l 


dr 


where 


dV(^±)i _ V2{^/aiRi + vj =F Vj) + T(±)i(i?i) + 7(±)i 

" 2(tT, ± 27/,) 


dT ^ 0-2 2viai 

Equations ( ^.41 ) and ( 2.42 ) gives the relation between uJoh{i) and uJohii — 1), for z > 2, as 

Wob(i) = f{i)iOoh{i - 1) , 


where 


/(*) = 


dV(_)i/dT 

dV{+),_ildT 


(2.41) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


For the first region, a direct calculation leads to 

<^ob(l) = — ^ (^Ri + +1^1^ = /(l)aJob(O) (2.46) 

where Wob(O) is the frequency of the light ray observed by an observer rest at the origin r = 0. Thus, using the above 
relations, we obtain the redshift of the light ray emitted from the z-th shell toward the observer rest at r = 0 in the 
form 


1 + z{i) 


^ob(») 
^ob (0) 


i 


n 

i=i 


(2.47) 


Our next task is to find the angular diameter-distance dA ■ The definition of dA is 

dA=j (2.48) 

where D is the physical length of the source perpendicular to the line of sight, and 6 is the observed angular size. 
Since the space we are considering is spherically symmetric, the diameter distance from the observer at the center to 
the z-th shell agrees with the circumferential radius Ri when the null geodesic intersects it; 


dA(z) = Ri 


We calculate dA-z relation in the dust-shell universe 


using the relations ( 2.39 ), (H) and ( ^.49| ). 


(2.49) 
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III. RESULTS AND DISCUSSION 


A. Setting the parameters of dust—shell model 

As mentioned in the previous section, the choice of the parameters and the initial hypersurface determines the 
behavior of a dust-shell universe. Since we are interested in cases which have a FL limit, we set the parameters so 
that they approach a FL model in the large N limit. 

For most cases described below, we set the mass distribution of shells as 

Pi = Pc (independent of i), (3.1) 

at the initial slice. This is not a unique but the simplest case which goes to a FL model in the large N limit. Using 
Pc, we define iLsheii and rn by 


/t2 _ -2 _ 

^shell — 'h — 


SttG 


(3.2) 


In terms of FL models, Lfsheii and rn may be regarded as the “Hubble constant” and “Hubble horizon radius” . 
However, we cannot say anything at this stage about what relation they have with the Hubble constant and horizon 
scale of FL models. 

For Xi, we put 


with a constant interval Aa; 


Ci('rinit) = fAx 




A^h 


where Ah is some positive integer. Here note that x^Vh = and the following relation holds 

2GM+(Ah) 


= 1 . 


XNh 


(3.3) 

(3.4) 

(3.5) 


Hence rn corresponds to the “mean Schwarzschild radius” of the An-th and (Ah + l)-t h reg ions. _ 

Be fore we proceed, we estimate the magnitude of the radiation-like term in equation ( 2.15 ). From equations (|2.13| ) 
and (2.8) with mi = 0, we find 


Gm2n-i = ]^(”- “ l)^(4n - l)rH, 

Gm2n = ^n^(4n - 3)rH, 

where n is a positive integer. From equations ( ^.16|) ( ^.14| ) and (^)^), we find msii) = mi+i — mi and obtain 

6n^ — 6n+ 1 


GTOs(2n - 1) =-^ 

Gms(2n) = 

A3 


-’’H, 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


Thus, when we consider a large A limit with fixing Xi = rn^i/Nn), the baryonic mass Gtos(j) = G(A^ ) is regarded 
as a small quantity, compared with GM+(f) and Gm^. That is, the radiation-like term is of order A^^ of the first 
term and can be neglected when Ah is large. 










B. Distance-redshift relation in “orthogonal” model 


As described in the previous section, the choice of Ei corresponds to choice of initial hypersurface. In a FL model, 
the simultaneous hypersurface is orthogonal to the trajectory of matter. Thus, we try choosing Ei so that the vector 
is orthogonal to the trajectory of each shell (we will refer to this choice as “model A”). From the condition 
— 1) = 0 at r = ri_i, we obtain 


A,; = 


n-i 


(Ji-lTi-l + 


(3.10) 


In Figure 2 we plot the angular diameter distance as a function of redshift of each shell for various A^h with a 
common i?sheU in the case of model A. The total number of shells A^t is taken to be 3 A^h- We emit a photon toward 
the center from each shell so that every photon reaches at the center simultaneously. We have obtained a surprising 
result that all the curves are quite similar even in the case when we put only two shells within the initial horizon 
radius; the deviation among the curves is at most about 10%. On the other hand, it can be also seen that the slope 
at higher redshifts becomes steeper as we decrease the number of shells. 

We compare the d^-z relation in the dust-shell universe obtained above with that of a FL m odel. We here adopt 
a spatially flat FL model, since the “curvature” term in the expansion law of a dust-shell (2.15) vanishes. Moreover, 
when A^h is large, the “radiation” term is negligible compared to the first term. Thus, in this paper, we focus on the 
cases with large A^h(> 10), and compare them with a spatially flat FL model filled only with non-relativistic matter. 
The cases with small 10) will be mentioned later. 

The Hubble equation in this flat FL model is 




(3.11) 


with pfl oc a~^. The relation between the redshift and the scale factor in FL models is 1 + z = oo/a where subscript 
“0” denotes a value when the observer receives the photon. Thus we can write the present Hubble parameter Hq in 
terms of the initial Hubble parameter Hi and redshift z; as 


H2 = iLf(l + zO-3 


(3.12) 


The dji^-z relation in the flat FL model is calculated once we fix Hi and zj, since the relation is determined solely by 
Ho- 

The dashed line in Figure 2 shows the c?a-z relation in the flat FL model with Hi = iLgheii- The redshift of the initial 
hypersurface is identified with the redshift of the outermost shell for the case A^t = 3A^h = 150 , i.e., Zi = z{i = 150). 
The deviation of the dust-shell universe from the flat FL model with = 150 amounts to about 2% around the 
maximum of the curve. We have confirmed that the deviation from the FL model becomes small as we increase A^h- 
For s maller A^h, however, the redshift of the outermost shell becomes small, resulting a larger Hq from the relation 


(I 


.12). We show in Figure 3 the dA-z relation in the dust-shell universe in the case of A^h = 10(solid line) and that 
in the flat FL model with Hi — iLsheii and Zi = z{i = 30) (long dashed line). In this case, the deviation amounts to 
about 10%. Moreover, the difference lies not only in the normalization of Hq, but also in the shape of the curve. To 
see this, we also plotted a FL curve (short dashed line) with a Hubble parameter changed by 10%: Hi = 0.9 x Hsheii- 
Note that the change in the Hubble parameter of the FL model only results in the change in the normalization of the 
curve. We see that the slope at higher redshifts is steeper in the dust-shell universe than in the FL model. It is also 
clear from Figure 2 that this tendency becomes strong as we decrease A^h- 

Thus, we can say that for large A^h, the flat FL model approximates the dust-shell universe quite well, but as we 
decrease Nh, this fit becomes poorer. 

In the next subsection we discuss the reason of this behavior by studying the behavior of averaged density, and try 
to reduce the deviation from FL models without increasing A^h- 


C. Behavior of Averaged density 

We usually regard a FL model as a large-scale “average” of a locally inhomogeneous universe. We will study the 
relation between the results obtained in the previous section and an “averaged” density of the dust-shell universe. 
We consider an averaged density around the observer at the center. The averaged density is defined by dividing the 
mass contained within some radius by the 3-volume on the hypersurface up to that radius. We will study the behavior 
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of the averaged density when the radius to take the average is gradually increased, and discuss its relation with the 

distance-redshift relation. _ 

First let us consider the 3-volume on the initial hypersurface defined in the previous section. From equation ( ^.30 ), 
the intrinsic metric of the initial hypersurface is given by 

- dr"^ + r'^dn'^. (3.13) 

Using the above line element, we obtain the spatial volume, Vol(*), of the *-th region (i > 1) on the initial slice in the 
form 


Vol(i) 


47r 


7* 2 


;_1 {1 + Ef)r - Hi 


dTT Y^(l + Ef)r - Hi 


dr 


r.5/2 


3(1+ U2) 


12il + Ef)^^ 8(1 + 


5TTHt^ , ( V+ E‘^)n - + V{I + E'^)r, \ 

2(1 + F;2)7/2 y i/(1 + Ef)ri_i - Hi + \/{^ + Ef)ri_i j 


(3.14) 

(3.15) 


For i = 1, Vol(l) is equal to 47rrf/3. 

Using this volume, we define the averaged density p{i) as follows. We sum the baryonic masses up to the (z — l)-th 
shell and half of the z-th shell, and divided the sum by the 3-volume inside the z-th shell; 

p{i)=\mS)/2+ ^ ms(j) W ^ Vol(fc) . (3.16) 

y j<2—1 J k<i 

It seems natural to take this sum of masses in averaging, since the motion of the z-th shell is approximately determined 
by sum of the gravitational mass inside the shell (which agrees with the sum of baryonic masses up to (z — l)-th shell) 
and half of its baryonic mass; t he numerator in the above definition is the same with M_|_(z) which appears in the 
first term in the expansion law (HI)- We plotted in Figure 4 the averaged density p{i) for the model described in the 
previous section. As expected, p{i) is almost constant near pc, which may explain the reason the deviation of d^-z 
relation in the dust-shell models from the FL curves is small. 

However, the averaged density slightly becomes smaller at the outside region than pc defined by the expansion law. 
The reason is explained as follows. The volume, U, between Vi and = Vi -\- Ar is expanded in terms of Ar as 

/*ri+Ar /*ri+Ar /Ti+Ar 

U = dzr / {F{r)}~^r'^dr = Att / {F{ri)}~^r^dr — Att / F(ri)'{F{ri)}~'^(r — ri)r^dr + ■ ■ ■ (3.17) 

Jn J n J n 

where F' = dF/dr. Substituting the expression of F, we can expand V in terms of 1 /Ah to find 

V = A7Tr^Ar{l + l^ + ---) . (3.18) 

Z iVH 

We see that the volume becomes larger than that of a homogeneous model. This effect is significant especially when 
z > Ah, i.e., beyond the horizon scale, which explains the behavior of the averaged density at large z in Figure 4. From 
this figure, one may think we can reduce the deviation from the FL model by adjusting the expansion law so that the 
averaged density is a constant value. We plotted the dA-z relation in this case (model B) and the corresponding FL 
curve for Ah = 10 in Figure 5 (the solid line and the long dashed line). We adjust pi in the expansion law iteratively 
so that the relation 


P(*) = (3.19) 

is satisfied. As a result, pi is not homogeneous; it increases as z increases (Figure 6). From Figure 5, we see that the 
difference between the dust-shell universe and the FL model remains the same as Figure 3, although the averaged 
density is indeed a constant. 
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D. Cases with homogeneous averaged density held 


We can take an initial hypersurface in which both pi and p{i) is homogeneous by requiring the relation 

Att 

Vol(*) = y(rf-rf_J . (3.20) 

We plot the d^-z relation of this model (model C) in Figure 7. Comparing with Figure 2, we see that the difference 
among the curves is reduced. Figure 8 illustrates the good agreement between the dA-z relations of the dust-shell 
model and a FL model for the case IVh = 10. 

Flere we check whether the assumption of constant interval Aa; = rn/iVu is essential or not. We try some patterns 
of the initial circumferential radius of the shells Xj, for instance, 

Xi+i = Xi + h X (3.21) 

where h and a are some positive constants and a 1. When a > 1(< 1), the distribution of shells becomes sparser 
(denser) for larger i. We have found that the results are unchanged. We also try another pattern 

X2n = 2n X Ax , (3.22) 

X2n+i = (2n + 1) X Ax + oAx , (3.23) 

with a 0. We plotted the result when Ah = 10, At = 30, and a = 0.05 (model D) in Figure 9^. We can see that 
the deviation from the FL model remains small. 

We can take another interesting choice; the expansion law is homogeneous (characterized by Pc), and the averaged 
density is also homogeneous, but each quantity differs; pc = ap (model E). This is realized by imposing a condition 
of the form 

A nr 

Vol(i) =ay(rf-rf_i) , (3.24) 

with a ^ 1. In Figure 10, we show the result in the case Ah = 10 and a = 0. 93 ( solid line). The dashed lines are 
curves in a FL model. The upper one uses the original “Hubble constant” (Eq.(^^, i.e., A, = Agheii, and the lower 
one uses the parameter changed by the same amount with the change in the volume, i.e., H? = x a~^. We 

see that the shape of the curve differs from a FL curve (compare with Figure 8). This indicates that the averaged 
density should agree with pc which determines the expansion law, in order for the dA—z relation to behave like that of 
a FL model. It should be also noted that the FL curve with the changed Hubble parameter is closer to the curve in 
dust-shell universe than the other. This may indicate that the “observed” Hubble parameter is closer to the averaged 
density p, rather than p^ which determines the expansion law. 


E. Summary of models and Conclusions 


Here we summarize our results. 

• Model A: The initial hypersurface is orthogonal to each trajectory of shells and the expansion law is homogeneous; 
Pi = Pc- The dA-z relation of the dust-shell universe shows deviation from the FL curve, especially when Ah 
is small (Figs.2,3). The averaged density p is not constant (Fig.4). 

• Model B: We adjust pi so that p = p^ is satisfied (Fig.6) on the same initial hypersurface as Model A. The 
deviation from the FL curve is not reduced (Fig.5). 

• Model C: We choose the initial hypersurface so that the relation p{i) = pi = Pc is satisfied. The dA-z relation 
of the dust-shell universe shows good agreement with the FL curve (Figs.7,8). 


“There is a maximum in a according to At when we impose the condition p{i) = pi = p^; if we increase At, the maximum 
value allowed for a becomes small. The value of a adopted here is about the maximum value for At = 3 Ah. 
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• Model D: The cho ice of the initial hypersurface is same with Model C, but the interval of the shells is not constant 
(Eqs. (|3.22| ),( 3.23 )). The deviation of cIa-z relation of the dust-shell universe from the FL curve remains small 
(Fig.9). 

• Model F: The averaged density p{i) and the parameter pi are homogeneous but p{i) ^ pi. The d^-z relation of 
the dust-shell universe shows mild deviation from the FL curve (Fig.10). 


From these results, we conclude that the d^-z relation in a dust-shell universe looks like a flat FL universe, when 
the expansion law resembles the flat FL model, and the behavior of averaged density field is scale-independent when 
we increase the scale of averaging, and the averaged density agrees with p^- This statement seems to be valid even 
in the cases with quite small number of the shells. However, the situation is not so simple. In small iVn cases, the 
radiation-like term in the expansion law cannot be neglected. One may expect that a FL model with radiation term 
gives a better fitting to those cases, but we have found this does not work. This implies that we cannot tell the effect 
of homogeneities just by studying the expansion law. We need more detailed study to this problem, which is left for 
our future work. 

We also note that in spatially flat cases the gravitational mass and the baryonic mass coincide; whether we use the 
baryonic mass or gravitational mass in defining the averaged density, the result is the same. In cases where those 
masses are different, we have to be careful in determining the averaged density when we try to construct a FL model 
which fits a dust-shell universe. In order to clarify which mass we should use to construct a fitting FL model, we 
have to study non-zero ki cases, which will also be done in our future paper. 


IV. SUMMARY 

We have studied the behavior of d\-z relation in a spherically symmetric dust-shell universe where the mass 
distribution is discrete. We have compared the relation of dust-shell universe with that of FL models, and discussed 
the relation with the behavior of averaged density. We have seen that the d^-z relation observed at the center agrees 
well with that of a flat FL model if the following conditions are satisfied: (i) the expansion law of the circumferential 
radius of the shells resembles the Hubble equation of a spatially flat FL model, (ii) the behavior of averaged density 
around the observer at the center is scale-independent as we increase the scale on which we take the average, and 
(iii) the averaged density agrees with the energy density of the FL model. We have also seen that the choice of initial 
hypersurface relates the expansion law to the averaged density. 

The effect of discreteness of mass distribution appears in the equation of motion of each dust-shell. This effect 
becomes smaller as we increase the number density of shell. We conclude the discreteness of matter distribution itself 
is of no significance in this model in discussing the observed quantities such as d^ and z, as long as the expansion law 
and the averaged density field is homogeneous in the sense described above. This supports the averaging hypothesis 
that a universe is described by a FL model if the universe is homogeneous when the density is averaged on some scale 
larger than the scale of the inhomogeneities. 

We need, however, further discussion for the cases when the number of shells is extremely small, and when the 
curvature term does not vanish. We also note that it will be interesting to study cases where the averaged density 
is inhomogeneous or non-radial null rays travel, to see cosmological lensing effects. These problems are left for our 
future work. 
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FIG. 2. Angular diameter distance-redshift relation in dust-shell universe for model A. Data points are connected by solid 
lines. The definition of the models is summarized at the end of section 3. The number of shells within the initial Hubble 
horizon A^h is 2, 5,10, and 50 from top to bottom around the maximum. The total number of shells At is taken to be 3Ah. 
We also mark the data by dots in the cases Ah = 2 and 5. We see that all the curves are quite similar; the deviation among 
the curves is at most about 10%. On the other hand, it can be also seen that the slope at higher redshifts becomes steeper as 
we decrease the number of shells. The dashed line shows the dA-z relation in the flat FL model with Hi = Asheii. The redshift 
of the initial hypersurface is identified with the redshift of the outermost shell for the case At = 150, i.e., Zi = z{i = 150). The 
deviation between the flat FL model and the dust-shell universe with At = 150 amounts to about 2% around the maximum 
of the curve. 
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FIG. 3. Angular diameter distance-redshift relation in the dust-shell universe for model A in the case of A^h = 10(solid line) 
and that in the flat FL model with Hi = /faheii and z\ = z{i = 30) (long dashed line). We see that the deviation among them 
amounts to about 10%. To see the difference in the shape of the curve, we also plotted a FL curve with a Hubble parameter 
changed by 10%: Hi = iLsheii x 0.9 (short dashed line). We see that the slope at higher redshifts is steeper in the dust-shell 
universe than in the FL model. 



i/Nh 

FIG. 4. Averaged density p normalized by pc in the cases A^h = 2, 5,10, and 50 for model A. We connect the data by solid 
lines. We also mark the data by dots in the cases A^h = 2 and 5. The curve which is smaller at large i corresponds to the curve 
with smaller A^h. We see that the averaged density p{i) is almost constant near pc, especially for small i, which may explain 
that the deviation of (Ia-z relation in the dust-shell models from the FL curves is small. One also notices that the averaged 
density becomes smaller at the outside region than pc. The interpretaion of this behavior is discussed in the text. 
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FIG. 5. Angular diameter distance-redshift relation in the dust-shell universe for model B in the case of A^h = 10(solid line) 
and that in the flat FL model with Hi = Hsheii and Zi = z{i = 30) (long dashed line). For comparison, we also showed the 
curve in the dust-shell universe for model A (short dashed line). We see that the difference in the normalization between the 
dust-shell universe and the FL model is reduced, but the difference in the shape of the curve is not reduced compared with 
Figure 3 . 



i/Nh 

FIG. 6. The averaged density p (dashed line) and the parameter pi (solid line) which determines the expansion law of the 
shells of model B. We see that pi increases as i increases, while the averaged density p is indeed a constant. 
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FIG. 7. Angular diameter distance-redshift relation in dust-shell universe for model C. The number of shells within the 
initial Hubble horizon A^h is 2,5,10, and 50. The total number of shells A^t is taken to be SA^h. Data points are connect by 
solid lines. We also mark the data by dots in the cases Ah = 2 and 5. Comparing with Figure 2, we see that the difference 
among the curves is reduced so that it is hard to distingush the curves for Ah = 5,10, and 50. 



Z 

FIG. 8. Angular diameter distance-redshift relation in the dust-shell universe for model C in the case of Ah = 10(solid line) 
and that in the flat FL model with Hi = Asheii and Zi = z{i = 30) (dashed line). This illustrates the good agreement between 
the cLa-z relations of the dust-shell model and of a FL model for the case Ah = 10. 


17 






z 

FIG. 9. Angular diameter distance-redshift relation in the dust-shell universe for model D in the case of A^h = 10(solid line 
and points) and that in the flat FL model with Hi = HsheW and Zi = z{i = 30) (dashed line). We can see that the deviation 
from the FL model remains small. 



Z 

FIG. 10. Angular diameter distance-redshift relation in the dust-shell universe for model E in the case of Ah = 10(solid line) 
and that in the flat FL model with Hi = iLsheii and Zi = z{i = 30) (short dashed line). The long dashed line is the FL curve 
with the Hubble parameter changed by the same amount with the change in the volume, i.e., H^ = iLgheii a~^. Gomparing 
with Figure 8, we see that the shape of the curve differs from a FL curve. This indicates that the averaged density should 
agree with pc which determines the expansion law, in order for the dA-z relation to behave like that of a FL model. Also note 
that the FL curve with the changed Hubble parameter (long dashed line) is closer to the curve in dust-shell universe than the 
other (short dashed line). The interpretaion of this behavior is discussed in the text. 
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